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standard model with U(2) flavor symmetry, in order to study physics potential of the present and 
QJ , future B factories. We evaluate CP asymmetries in various B decay modes, AmB^, Am^^, and 



We investigate quark flavor signals in three different supersymmetric models, the minimal su- 
pergravity, the SU(5) SUSY GUT with right handed neutrinos, and the minimal supersymmetric 



ek- The allowed regions of the CP asymmetry in B J/ij^Ks and AmB^/ AmB^ are different 



^ . for the three models so that precise determinations of these observables in near future experi- 
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5t , ments are useful to distinguish the three models. We also investigate possible deviations from 

the standard model predictions of CP asymmetries in other B decay modes. In particular, a 
large deviation is possible for the U(2) model. The consistency check of the unitarity triangle 
including B vrvr, pir, D^*^K^*\D^*K,D*p, and so on, at future high luminosity e^e B facto- 
ries and hadronic B experiments is therefore important to distinguish flavor structures of different 
supersymmetric models. 
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I. INTRODUCTION 



Recent development in B physics is remarkable. Both the Belle experiment at KEK and 
the BaBar experiment at SLAC observed large CP violation in B ^ J /if) Ks and related 
modes. These observations are the first discovery of CP violation out of the kaon system[0, 
The results are consistent with the Kobayashi-Maskawa mechanism of CP violation in the 
three-generation Standard Model (SM). In the coming years, we expect much improvements 
in measurements of CP violation and rare decay processes at the asymmetric B factories. 
In addition, the magnitude of the B^ — Bg mixing will be determined at the Tevatron 
experiments!^]. It will be clear in a few years whether or not the Cabibbo-Kobayashi- 
Maskawa (CKM) matrix is the main source of the flavor mixing and the CP violation in the 
quark sector. 

In future, B physics is expected to play an even more important role in precisely deter- 
mining the flavor structure of the SM and searching for possible new physics effects beyond 
the SM. LHC-B|^ and BTeV0 experiments are planned to provide very precise information 
on the angles of the unitarity triangle from B^ and Bg decays at hadron machines. As for 
e~^e~ asymmetric colliders, both KEK and SLAC are considering to increase luminosity by 
1 — 2 orders of magnitude by the time that these hadron experiments will be carried out[0, 
With luminosity of 10^^ — 10^^cm~^s~^, the super B factories will provide us with 10^ — 10^° 
BB pairs in a year. Then, we shall have good opportunities to explore new physics from 
observations of CP asymmetries and rare B decay processes. 

Among various candidates of new physics beyond the SM, supersymmetry (SUSY) is the 
most interesting one. Although the main motivation for introducing SUSY is to solve the 
hierarchy problem, namely to give a justification to the electroweak scale, which is much 
smaller than the Planck scale, flavor physics can provide important information on SUSY 
models. In SUSY models, mass matrices of SUSY partners of the quarks and the leptons 
are new sources of the flavor mixing. Since these mass matrices are determined from SUSY 
breaking terms in the Lagrangian, their flavor structures reflect the SUSY breaking mecha- 
nism and interactions present between the scale of the SUSY breaking and the electroweak 
scale. Future B physics is therefore very important to discriminate various SUSY breaking 
scenarios. It can play a more important role if LHC experiments discover SUSY particles, 
in which case we can make more precise predictions for flavor signals based on a particular 
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scenario. 

In this paper, we investigate SUSY effects on B physics based on three different SUSY 
models, namely (1) the minimal supergravity model (mSUGRA), (2) the SU(5) SUSY GUT 
with right-handed neutrinos, and (3) the SUSY model with U(2) ffavor symmetry]^, |TU|. 
We focus on the new physics search through the consistency test of the unitarity triangle. 
We address two questions. First, we ask whether these models can be distinguished from 
the SM in near future by measuring the CP asymmetry of B —>■ J/ipKs and the Bg — Bg 
mixing. Second, we consider impacts of other angle measurements of the unitarity triangle 
in the era of LHC-B/BTeV and an e^e~ super B factory when the above two observables 
are precisely determined. We analyze the three models in the same fashion, so that we can 
clearly show the potential of future B physics. There are already many analyses in the 
literature for each of three models[^ |10|, |ll], |12|, |13|, but here we make a systematic 



treatment to show that various measurements in B physics are in fact useful to distinguish 
different SUSY models. We show that the allowed region of AtubJ Auib^ is quite different 
from the SM prediction for the SU(5) SUSY GUT with right-handed neutrinos and the U(2) 
model, whereas the deviation is small for the mSUGRA. Furthermore, the GUT and the 
U(2) model can be distinguished when we measure CP asymmetries of decay modes such as 
B -> Tin, pTT, D^*^K^*\ D^*h, and D*p. 

This paper is organized as follows. The three models are introduced in Sec. II. The 
Bd — Bfi mixing, the Bg — Bg mixing, the CP violating parameter in the — mixing 



{sk)-, and CP violations in various B decays are discussed in Sec. |Tl[ The numerical results 



of these observables are presented in Sec. |V[ Our conclusion is given in Sec. 0. 

II. MODELS 

A. The Minimal Supersymmetric Standard Model 

The minimal supersymmetric standard model (MSSM) is an SU(3)c x SU(2)l x U(1)y su- 
persymmetric gauge theory with the SUSY being softly broken. The MSSM matter contents 
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are the following cliiral superfields: 

Q,(3, 2, i) , I7,,(3, 1, -h , A(3, 1, \) , 

O 6 6 

2, , 1, 1) , 

^i(l,2, -^), i^2(l,2, i), (1) 

where the gauge quantum numbers are shown in parentheses and i = 1, 2, 3 is a genera- 
tion index. Assuming R-parity invariance and renormalizability, we can write the MSSM 
superpotential as 

Wmssm = fj^D,Q,Hi + fjU,Q,H2 + fj^E^LjE, + ^iH^H^ . (2) 
The soft SUSY breaking is described by the following Lagrangian: 

- Aoft = + {ml)/d% + {ml\'^% 

+ (m|)Ve,?^^- + K)/Ttf, 
+A\h\hi + Alh\h2 - {Bfihih2 + H.c.) 
+ {^A^diqjhi + A'^Uiqjh2 + Afeiljhi + H.c, 

+ '^BB+'^WW+'^gg, (3) 

where g,, Ui, di, li, Cj, hi, and /12 are the corresponding scalar components of the chiral 
superfields, and B, W, and 5^ denote U(1)y, SU(2)l, and SU(3)c gauge fermions, respectively. 



B. Flavor Structure of the Soft Breaking Terms 

Although the Yukawa couplings are the only source of the flavor mixing in the SM, the 
mass terms and the trilinear scalar coupling terms (A-terms) of squarks and sleptons in 
Eq. may induce additional flavor mixings in the MSSM. 

The Yukawa couplings /'s in Eq. (§) are constrained to reproduce the known quark 
and lepton masses and the CKM matrix. On the other hand, the soft breaking terms, 
their flavor structures in particular, are rather unconstrained at first sight, apart from the 
naturalness condition that they should be within the TeV scale. As is well-known, however, 
unless some specific structure is assumed in the soft breaking terms, the sfermion-exchanging 
contributions to flavor changing neutral current (FCNC) processes such as the — 
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mixing and the /i — 67 decay are too large to satisfy the experimental limits, if the squark 
and slepton masses are below a few TeV. There are several ways to avoid this problem. 

One is to assume a SUSY breaking (and its mediation) mechanism in which the uni- 
versality of the soft breaking terms are ensured. In other words, mass degeneracy for the 
sfermions with the same electric charge and chirality, and proportionality of the A-terms to 
the Yukawa coupling constants are required to suppress FCNC processes. Phenomenology 
of models with the universality further depends on the energy scale where the SUSY break- 
ing is generated, because the universality in the sfermion sector is vitiated due to radiative 
corrections induced by the Yukawa couplings. 

It is convenient to use renormalization group (RG) equations in order to trace these 
radiative corrections. The universality in the soft breaking terms is imposed on the boundary 
conditions of the RG equations at the energy scale of the SUSY breaking. If the energy scale 
of the SUSY breaking is close to the electroweak scale, the RG evolution is tiny and the 
universality is practically maintained. The gauge mediated SUSY breaking model [0 is of 
this kind, and its flavor phenomenology is quite similar to that of the SM. 

If the energy scale of the SUSY breaking is far above the electroweak scale, the RG 
evolution is sizable and the universality is lost. Though not exactly universal, flavor physics 
is still under control in this kind of models in the sense that the origin of the flavor mixing 
only resides in the Yukawa couplings. The flavor mixing in the squark sector is determined 
by the quark masses and the CKM matrix. On the other hand, the flavor mixing in the 
slepton sector is ruled by the lepton couplings in the superpotential including Majorana 
mass terms of right-handed neutrinos if exist. The mSUGRA discussed in Sec. |II Q is a 
model in this class. 

Embedded in a GUT, the above situation is modified if the energy scale of the SUSY 
breaking is higher than the GUT scale. Since the GUT interactions obscure the distinction 
between the quark flavors and the lepton flavors, the lepton flavor mixing in the Yukawa 
couplings affects the squark sector. We shall examine an SU(5) SUSY GUT with right- 
handed neutrinos among models of this kind in Sec. [II D| . 

Another way to suppress the FCNC processes is to rely on a flavor (or horizontal) sym- 
metry. It is obvious that the flavor symmetry, whatever it is, should be broken because the 
Yukawa couplings have no such symmetry. The symmetry breaking must be taken so that 
the observed quark and lepton masses and their mixings are reproduced. Even though this 



6 



constraint is imposed, there are several choices for the flavor symmetry and its breaking 
pattern. The flavor phenomenology heavily depends on them. In Sec. [11 E| , we shall consider 
a model with U(2) flavor symmetry among the possibilities. 



C. The Minimal Supergravity Model 

The mSUGRA consists of the observable sector, i.e. the MSSM, and a hidden sector. 
These two sectors are only interconnected by the gravitation. The SUSY is assumed to 
be spontaneously broken in the hidden sector, and the soft breaking terms in Eq. (0) are 
induced through the gravitational interaction in the following manner: 

Al = Al = ml, 

A% = moAofi , A'l = moAofJ , A'l = m^A^fi , 

Ml = M2 = M3 = Ml/2 , (4) 

where we assume the GUT relation among the gaugino masses. The above relations are 
applied at the energy scale where the soft breaking terms are induced by the gravitational 
interaction. We identify this scale with the GUT scale {Mx — 2 x lO^^GeV) for simplicity. 

The soft breaking terms at the electroweak scale are determined by solving RG equations 
with the initial conditions defined in Eq. (|^). 



D. The SU(5) SUSY GUT with Right-handed Neutrinos 

The measurements of the three gauge coupling constants at LEP, SLC, and other ex- 
periments support the idea of the supersymmetric grand unification. Furthermore, there 



is clear evidence of neutrino oscillations in the atmospheric [0 and the solar ||18[ neutrino 
experiments. Guided by these experimental results, SU(5) SUSY GUT with right-handed 
neutrino has been studied. In particular, the relationship between quark flavor signals and 
the neutrino oscillations has been investigated in Ref. |jl2|, |13|, |I4|. A large flavor mixing 



in the neutrino sector can include a squark mixing in the right-handed down-type squark 
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sector. In Ref. ||T2[, the quark flavor signals are studied for various neutrino oscillation sce- 
narios. In Ref. |T3[, effects of CP violating phases in the GUT Yukawa coupling constants 
are taken into account. It is shown in these papers that large contributions to Ek and the 
fi ej decay can arise from the new source of flavor mixing in the neutrino sector. These 
analyses are extended to the case of a GUT model with realistic fermion mass matrices in 
Ref. []T1|. Here we follow the analysis of Ref. []T^| and we give a brief description of the 



model. 

The Yukawa couplings and Majorana masses of the right-handed neutrinos in the model 
are described by the following superpotential: 

o 

+i\NrWFj)aH'^ + ^{M^ymjv^, (5) 

where i and j are generation indices, while a, b, c, d, and e are SU(5) indices, eabcde denotes 
the totally antisymmetric tensor of the SU(5). Tj, Fi, and Ni are 10, 5, and 1 representations 
of the SU(5) gauge group, respectively. Tj contains Qi, Ui, and Ei in Eq. (|I|), and Fi includes 
Di and Lj. H and H are Higgs superfields in 5 and 5 representations, respectively. H 
consists of Hc{S, 1, — |) and H2, and H does Hc{S, 1, |) and Hi. {XuY^, (Ad)*-', and {XnY^ 
are the Yukawa coupling matrices, and {M^y^ is the Majorana mass matrix. In addition to 
the above superpotential, we also need a superpotential for Higgs superfields, Wh{H, H, E), 
where S"^ is a 24 representation of the SU(5) group. It is assumed to develop a vacuum 
expectation value (VEV) as (S"^) = diag(|, |, |, — |, —^)vg at the GUT scale and breaks 
the SU(5) symmetry to SU(3)c x SU(2)l x U(1)y. 

The supermultiplets whose masses are of order of the GUT scale such as He and He can 
be integrated out below the GUT scale. Then, the effective theory below the GUT scale is 
the MSSM with the right-handed neutrino supermultiplets, and its superpotential is given 
as 

Wmssm.^ = Wmssm + {fNfNrL,H2 + ]^{MNfNj^j, (6) 

where the Yukawa coupling matrices are related to those in Eq. as {fuY'^ = (Ac/)*"', 
{.foY^ = {IeY^ = (Ad)*"', and {/nY'' = {^nY'' in the leading order approximation. 

In the energy region lower than the Majorana mass scale (= M^), the singlet supermul- 
tiplets are integrated out, and the resulting superpotential is the sum of Wmssm in Eq. (^ 



and the following higher dimensional term: 

AW. = -\{K^f{UH,){L,H,), = (fj;y\±-),,(f^y^. (7) 

This term yields the neutrino masses below the electroweak scale as 

{m^r = {K,f{h2f- (8) 
The above neutrino mass matrix is related to the observable neutrino mass eigenvalues and 



the Maki-Nakagawa-Sakata (MNS) mixing matrix |[T9[] as 



- (^MNs)^^/(nlNs).' (9) 



in the basis in which the charged lepton mass matrix is diagonal. 
As is mentioned above, the superpotential in Eq. (H) predicts that 

UeT = UdT (10) 

at the GUT scale. It is well-known that the mass ratios of the down-type quarks to the 
charged leptons in the first and the second generations can not be explained by this relation 
although it reasonably works for the third generation. This, however, is not a fatal flaw of 
SU(5) GUTs because there are several ways to overcome this shortcoming. For example, 
higher dimensional operators with S^^, may contribute differently to the Yukawa coupling 
matrices of the down-type quarks and the charged leptons. 



In Ref. [0], quark FCNC processes, lepton flavor violation processes, and the muon 
anomalous magnetic moment were calculated in this model. To account for quark and 
lepton mass ratios, the following higher dimensional operator was introduced: 

AWsu{5).. = ^^maT.\{T,rH,, (11) 

where Mp is the Planck mass (Mp ~ 2 x 10^^ GeV). Consequently, Eq. ( p!0[ ) is modified to 

UeT = UdY - hi^oT. (12) 



where ^ = vajMp ~ 0.01. Taking the Majorana mass matrix proportional to the unit 
matrix ((Mat)*-' = MpS'^^) for simplicity, Baek et al. [T^ showed that the flavor mixings 



of the squark and slepton sectors were determined by the CKM matrix, the MNS matrix, 
and two additional mixing matrices related to the down-type quark and the charged lepton 



Yukawa coupling constants. As long as we take the large mixing MSW solution for the solar 
neutrino anomaly, the SUSY effect becomes large for ek and B(/i ^ 67). In this paper we 
follow Ref. but consider the special case that the two additional mixing matrices are 



equal to the unit matrix, because the general features mentioned above do not change by 
this simplification. 

In order to calculate the FCNC processes we need to specify the soft breaking terms. The 
SU(5) invariant soft breaking terms are written as 

o 



1 



{F )j:\{T^f'H, + H.C. 



+^M5G'5G5, (13) 

where T\ F*, and are the scalar components of T*, F*, and iV*, respectively; H and 
H stand for the corresponding scalar components of the superfields denoted by the same 
symbols; and G5 represents the SU(5) gaugino. We assume that the soft breaking terms are 
universally generated at the Planck scale, i.e. 

(m^)/ = (m|)/ = (m^)/ = mlSl 
{Xy^ = moAo(A)^^ (A = At;,Az5,Ajv), 

Ms = Ml/2. (14) 

These equations serve as a set of boundary conditions of RG equations at the Planck scale. 

We solve the RG equations of the SU(5) SUSY GUT from the Planck scale to the GUT 
scale, then those of MSSM with right-handed neutrinos between the GUT scale and M^. 
Finally, the squark and slepton mass matrices are obtained by the RG equations of the 
MSSM below Mr- 
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E. A Model with U(2) Flavor Symmetry 



It is possible that the family structure of the quarks and the leptons is explained by some 
flavor symmetry. Although U(3) is a natural candidate of the flavor symmetry, it is badly 
broken by the top Yukawa coupling. It is therefore legitimate to choose a U(2) symmetric 
model in order to study the flavor problem in the MSSM. 

In this framework, the quark and lepton supermultiplets in the first and the second 
generations transform as doublets under the U(2) flavor symmetry. Each of these doublets 
carries a positive unit charge of the U(l) subgroup. The quark and lepton supermultiplets 
in the third generation and the Higgs supermultiplets are totally singlet under the U(2). 
In addition to the ordinary matter fields, we introduce the following superfields: a doublet 
$*(—!), a symmetric tensor 5'*-'(— 2), and an antisymmetric tensor A^^{—2), where i and j 



run from 1 to 2, and the numbers in the parentheses represent the U(l) charges [p!0 . 

The U(2) invariant superpotential relevant to the quark Yukawa couplings is given as 
follows: 

>Vu(2) = yu (U3Q3H2 + ^W,QsH2 + ■^Us<l>'Q,H, 

+^SmQ,H, + ^A^m^Q.H,^ (15) 
+yD (dsQsH, + ^<i>W,QsH, + ^Ds<l>'Q^H, 

\ IVlp IVlp 

where Mp is the scale of the flavor symmetry, and yg, ag, Bq, cq, and dg {Q = U,D) 
are dimensionless coupling constants. In Eq. (|15|), we neglected dimension five and higher 
dimensional operators in the superpotential. Absolute values of the above dimensionless 
coupling constants except for ypi are supposed to be of 0(1). 
The successful breaking pattern of the U(2) symmetry is that 

U(2) — > U(l) — > l(no symmetry), (16) 

where the first breaking is induced by VEV's of and and a VEV of A'-' brings about 
the second one. The VEV's are given as 

Mp Mp Mp ^ ^ 
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where e$ and e' are taken to be real without loss of generality. Note that {S'^^) is chosen 
so that it leaves a U(l) unbroken. Hierarchical relations among the VEVs that e' <^ 6$ ~ 
1 65 1 ^ 1 are assumed in order to reproduce the quark masses and the quark mixing angles. 
With the above VEVs, we obtain the following quark Yukawa couplings: 



VQ 



\ 



Q = U,D, 



(18) 



age' 
-age' dge bge 
cqe 1 J 

where we use e = e<i, = €5, which is valid providing appropriate redefinitions of the coupling 



v 



constants in Eq. (|l^). Eq. (|18D is applied at the GUT scale where we assume that the 
symmetry breaking sequence in Eq. (16) occurs. The parameters in Eq. ([T8|) are determined 
so that the known quark masses and mixing are reproduced taking the RG evolution into 
account. 

The U(2) symmetry constrains the soft breaking terms as well as the supersymmetric 
terms. The U(2) invariant soft breaking terms relevant to squark masses are 



-c. 



E 



m, 



/2 



-,SS 



(19) 



where a/'s are dimensionless coupling constants of 0(1), and shown are the terms that yield 
squark masses of O(e^) or larger when the flavor symmetry breaking takes place. The squark 
mass matrices stemming from Eq. (^) are parameterized as: 



m 



X 



m, 



X2 



( 



1 

1 




Y 9 Y 

Toot rt-^t 



22 

'23 ^ 



23"= 
33 



X = Q, f/, D, 



(20) 



where r"^'s are constant parameters of 0(1). 

As for the A-terms, it turns out that they have the same hierarchical structure as the 
Yukawa couplings in Eq. (|T8|) : 



( a^e' ^ 

-age' dqt bqe 
cqe 1 J 



Q = U,D. 



(21) 
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In general, though being of 0(1), ag, bg, cq, and dq take different values from the corre- 
sponding parameters in Eq. (|1^). Therefore, we expect no exact universality of the A-terms 
in this model. 

The soft SUSY breaking terms at the electroweak scale are given by solving the RG 
equations of the MSSM with the boundary conditions in Eq. (^) and Eq. ( pi]) at the GUT 
scale. 

III. OBSERVABLES 

The observables considered in the following are the CP violation parameter ek in the 
mixing, Bd — Bd and Bs — Bg mass splittings Am^^ and Am^^, respectively and 
CP asymmetries in various B decay modes. 

The Bd — Bd, Bg — Bg, and — mixings are described by the effective Lagrangian 
of the following form: 

>C = CLLirLrQLa)ifLl,QLp) + CRniqRl'QRa){fRl,QRp) 
+Cili(tRQLe.)i4QLp) + C^Lli<lRQLp){4QLa) 

+C^Rl{qlQRo.MQRp) + C^RiqlQRmlQRc.), (22) 

where (g, Q) = {d, b), (s, b) and {d, s) for the Bd — Bd, Bg — Bg, and — mixings, 
respectively. The suffices a and (3 are color indices. The Wilson coefficients C's are obtained 
by calculating box diagrams. See Ref. for explicit formulas of the coefficients. The 
mixing matrix elements Mi2{Bd), Mi2{Bs), and Mi2{K) are given as 



Mi2(P) = 




where P = Bd, B^, K^. 

In the SM, the flavor changes only occur in the interaction with the left-handed quarks, so 
that Mi2 is dominated by the Cll term for all the three mesons. The situation is the same in 
the mSUGRA, since the flavor mixing in the squark sector is induced by the running effect 
and hence takes place only in the left-handed squark sector. In the SU(5) SUSY GUT with 
right-handed neutrinos and the U(2) model, on the other hand, there are sources of squark 
flavor mixing other than the CKM matrix. In the SU(5) SUSY GUT with right-handed 
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neutrinos, flavor mixing in the right-handed down-type squark sector is induced due to the 
Yukawa couphng matrix of the neutrinos through the running between the GUT and the 
Planck scales. In the U(2) model, the squark mass matrices contain more free parameters. 
Consequently, flavor mixing is possible in both the left-handed and the right-handed squark 
sectors and the mixing matrices can be different from the CKM matrix. In the latter two 
models, all the Wilson coefficients in Eq. (|22|) are relevant. 

We parameterize the matrix elements of the operators in Eq. (^) as 

{K'\{dirsLa){4l,SL^)\K') = (24a) 



{K\dls,^){dlsnp)\Kl = rny^B^^^^'', (24b) 

D \TrLg -\- 111(1 J 

{KXdlsno){dlsnm') = f^^T^V ^/^^^^ (24d) 



{Krdlsn,){dlsn.)\K') = 1 f^^V m^/^E^^ (24e) 

iz \ms + rud / 

where Bk, B^^^^''^\ and B^^^^''^^ are bag parameters of 0(1), which have been calculated by 
the lattice QCD method It can be seen that the matrix elements of the scalar operators 
are enhanced by a factor ~ {mK/fris)'^ for the — mixing. For B^ — B^ mixing cases, 
the corresponding factor is ~ {mB/mb)'^ so that the enhancement is less significant. In the 
SU(5) SUSY GUT with right-handed neutrinos and the U(2) model, the Clu and/or Cjiji 
terms can significantly contribute to Mi2{K) because of this enhancement in the matrix 
elements. 

We can express Ek, Am^^, and Ams^ in terms of M12 as 

e--/^Im Mi2(i^) ^25) 
\/2AmK 

AniB, = 2\M,2{Bd)\, (26) 

AniB, = 2\Mu{Bs)\. (27) 

The CP asymmetry in i? — > J/ipKs, Aqp{B — > J/tpKs), is defined as 
r{Bd{t) ^ J/4jKs) - T{Bd{t) ^ J/4jKs) 



TiBdit) ^ J/^PKs) + TiBdit) ^ J/^jjKs) 
This asymmetry is given by the phase of Mi2{Bd) as 



-A'^^iB J/^ Ks) sin Am^.t. (28) 



A5^?(5 ^ J/^ Ks) = sin ct^M , (29) 
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FIG. 1: The unitarity triangle. 

where 0m is defined as e"^'^' = Mui^Ba) /\Mi2{Bd)\. In the present analysis we assume that 
tree-level diagram dominates the Ba{Ba) J/ipKs decay so that no new phase appears in 
the decay amplitude^ . Experimentally, sin ipM can be determined by combining decay modes 
with the h — > ccs transition such as Bd J/ipKs, B^ —>■ J/ip Kl, and Bd — > ip' Ks- 

In order to constrain new physics from the consistency check on the closure of the unitarity 
triangle, depicted in FIG. |l], it is important to measure the angles other than (pM- There 
are several theoretically clean ways to determine these angles. In the SM, 20i is given 



by (pM, and sin 202 is obtained from the isospin analysis of B nn\21\ and the time 



dependent Dalitz analysis of -B ^ p7r||22[. B modes provide us with the angle 

03|23[. B D^*^7i^^^ and B D*p\25\ with the angular analysis give us information on 



sin(20i + 03). 

If we consider effects of new physics, these measurements can be interpreted as constraints 
on the phases of the Bd — Bd mixing and the decay amplitudes. When we neglect the new 
phases in the decay amplitudes, B ^ rni, pn, B D'-*^K^*\ and B D^*\, D*p provide 
us with sin(0A/ + 203), 03, and sin(0A/ + 03), respectively, where 03 is the weak phase of the 
h ^ u transition amplitude in the standard phase convention given in APPENDIX In 
the following analysis, we assumed that 0a/ is determined from the B —>■ J/ipKs mode and 
related modes, and we study impacts of the 03 measurement on new physics search. 



^ In the U(2) model, in particular, there might be sizable contributions to the decay amplitudes with new 

CP phases. In such a case a direct CP asymmetry may be observed. 
^ This approximation is valid for the three models under consideration, at least for B 

D'-*'>K^*\ D'^*^'k, D* p. New phases could be important for the decay amplitudes of _B ^ tttt, pir. Even in 

such a case, we could obtain information about new phases by measuring CP asymmetries of the various 

modes listed above. 
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IV. NUMERICAL ANALYSIS 



A. Parameters in the minimal supergravity model 

In our calculation, the masses and the mixing matrices in the quark and lepton sectors 
are treated as input parameters which determine the Yukawa coupling matrices. 

The CKM matrix elements Vus, Vcb, and \Vub\ are determined in experiments indepen- 
dently of new physics contributions because these are extracted from tree-level processes. 
We fix Vus and Vcb in the following calculations as Vus = 0.2196 and Vcb = 0.04, and vary 
\Vub\ within a range \Vub/Vcb\ = 0.09 ± 0.01. Although the current error of \Vub\ is estimated 
to be larger than this value, we expect theoretical and experimental improvements in near 
future. We vary the CP violating phase, 03, within ±180° because it is not constrained by 
the tree-level processes independently of new physics contributions. For the quark masses, 
we take m^^" = 175 GeV, m^"^" = 4.8 GeV, mP°i^ = 1.4 GeV, and mf^{2 GeV) = 120 MeV. 

As for the SUSY parameters, we assume M1/2, Aq, and fi are real parameters in order to 
avoid too large SUSY contributions to the electric dipole moments of the neutron and the 
electron. We vary these parameters within the ranges < rriQ < 3 TeV, < M1/2 < 1 TeV, 
and —5 < Aq < +5 at the GUT scale. Both signs of fi are considered. We take the ratio of 
two VEV's tan/3 = (/i2)/ (^1) = 20 for most of our analysis and comment on other cases. 



B. Parameters in the SU(5) SUSY GUT with right-handed neutrinos 



In the SU(5) SUSY GUT with right-handed neutrinos, we need to specify the parameters 
in the neutrino sector in addition to the quark Yukawa coupling constants given in Sec. IV A. 



We take the neutrino masses as ml^ — ml^ = 2.4 x 10 ^ eV^, m^, — ml^ = 4.2 x 10 ^ eV^, 



1^3 

and rriiy^ ~ 0, and the MNS matrix as 



1^2 



( 



V. 



MNS 



CsolCls SsolCia Si3 

~"5solCatm CsolSatm'^lS CgQiCatm '^solSatm'^lS ^aXrcfiVi 
'^sol'^atm CsolCatm'^lS Csol-^atm '^solCatm'^lS ^atmClS 

,2 



(30) 



(cj = cos^j. Si = sin 6^4) with sin^26'atm = 1, sin^ 26'soi = 0.655, and sin^26'i3 = 0.015. These 
mass differences and mixing angles are consistent with the solar and atmospheric neutrino 
oscillations. The value of sin^ 26'i3 is constrained by reactor experiments p6[, and the above 
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value is take as an illustration. 

In addition, we assume that the mass matrix of the right-handed neutrino in Eq. (|]) is 
proportional to the unit matrix, and we take = 4 x 10^^ GeV. Complex phases in the 
MNS matrix and the right-handed neutrino mass matrix are neglected. 

The soft SUSY breaking parameters in this model are assumed to be universal at the 
Planck scale, and the running effect between the Planck and the GUT scales is taken into 
account. We scan the same ranges for mo, M1/2 and as those in the mSUGRA case. 



Parameters in the U(2) model 



In the U(2) model, we take the symmetry breaking parameters e and e' as e = 0.04 and 
e' = 0.008, and the other parameters in the quark Yukawa coupling matrices in Eq. (|I8]) are 
determined so that the CKM matrix and the quark masses given in Sec. [IV A| are reproduced. 
Note that the texture of the Yukawa coupling matrices in Eq. (^) predicts the following 
relation among quark masses and CKM matrix elements: 

(31) 
(32) 



ruu _ 




2 


r 


Vub 








rric 






rrid _ 


Vtd 


2 


1 + 


Vtd 






rris 


Vts 




Vts 



In the numerical calculation we adjust m„ and md to satisfy these relations. 
There are many free parameters in the SUSY breaking sector as shown in Eq. 



OD and 



Eq. ([21|). In order to reduce the number of free parameters for numerical calculations, we 
assume that 



m, 



Q2 



m, 



U2 



m^2 = ml, 



' ij ~ ' ij — ' *J ' 



(zj) = (22),(23),(33). 



(33) 
(34) 



We vary these parameters within the ranges < mo < 3 TeV, — 1 < < +1, < < 4, 
|r23| < 4 and —180° < argr23 < 180°. The boundary conditions for the A parameters and 
the slepton mass matrices are assumed to be the same as the mSUGRA case to simplify the 
numerical analysis. We think that the above assumptions on the soft breaking terms are 
sufficient for our purpose of comparing new physics effects related to the 5° — and the 
— mixings in the three models. 
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D. Experimental constraints 



In order to obtain allowed parameter regions, we impose the following experimental con- 
straints: 

• Lower limits on the masses of SUSY particles and the Higgs bosons given by the direct 



search in collider experiments ||27|] 



Branching ratio of the 6 — 57 decay: 2 x 10 ^ < B(6 ^ 57) < 4.5 x 10 . 



Branching ratio of the fj, ^ decay for the SUSY GUT case: B(/i ^67) < 1.2 x 
10-"|2|. 



Measured values of Ek and AmB^[30|, and the lower bound of Am^JSl 



• CP asymmetry in the B —>■ J/ip Kg decay and related modes observed in the B factory 
experiments |]1|, 0. 

Although the values of sk and ^frtB^ are precisely measured in experiments, there are 
theoretical uncertainties in the evaluation of the matrix elements for /S.S{B) = 2 operators. 
In order to take these theoretical uncertainties into account, we calculate Ek and ^fnB^ with 
bag parameters and /^^ ^ in Table | and allow parameter sets if the calculated values of Ek 
and AniBj lie within the ranges 

= e*''/^(2.28 X 10"^) X (1 ± 0.15), (35) 
ArriB^ = 0.479 ps"^ x (1 ± 0.2)1 (36) 

For ArriB,, we impose a constraint on the ratio to Auib^ as AtubJ ArriB^ > 27 since a large 
portion of the theoretical uncertainties is expected to cancel by taking the ratio. For the 
CP asymmetry, we use A^p{B J/ipKs) > 0.5. 



E. Numerical results 

At first we discuss qualitative features of the SUSY contributions to the B^ — B^ mixing, 
the Bs — Bg mixing and Ek for each model. 

In the mSUGRA, it is well-known that the main SUSY contributions to Mi2(-Bd), 
^12 (-Bs), and Mi2{K) come from the box diagrams with the charginos and the up-type 
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fx 


Bk 




„LiJ{l) 


„LR(2) 




f.RR{2) 


159.8 MeV 


0.69 




1.03 


0.73 


0.65 


1.05 


fBa 


fSs/fBd 


Bb 


BdtBs 


Bd^Bs 






210 MeV 


1.17 


0.8 


0.8 


0.8 


0.8 


1.19 



TABLE I: Decay constants and bag parameters for the B^ — and the — mixing ma- 
trix elements used in the numerical calculation |2C]. are given by B^^^'^^ = 5B^^^^^ — 



^1 



4( — r^) ^Br with q = d,s, which is valid in the static limit. 

squarks and that the flavor mixing in the chargino vertex is determined by the CKM matrix. 
Consequently, the SUSY contributions are approximately proportional to the CKM matrix 
elements {V^^Vaf, {V*Vtby, and (V;^^^,)^ for Mi2(5,), M.^iB,), and M.^iK), respectively, 
and the ratios to the corresponding SM contributions are common: 

(AmeJsusY _ (Am^JsusY _ (gx)susY 
(Am^JsM (AmijJsM {£k)sm 

For the other two models, this proportionality is violated due to the squark flavor mixing 
induced by sources other than the CKM matrix. 

In the SU(5) SUSY GUT with right-handed neutrinos, the flavor mixing in the right- 
handed down-type squark sector is related to the flavor mixing in the left-handed slepton 
sector, and hence the constraints from the lepton flavor violating processes such as — > e 7 is 
important. As shown in Ref. [0], B(yU -^67) exceeds the present experimental upper limit 
1.2 X 10^^^ in the parameter region where the SUSY contributions to the 5° — mixings 
become larger than ~ 10%^. Therefore ArriBd, Am^^ and Aqp{B J/ipKs) are almost 
the same as the SM values for a given CKM matrix. On the other hand, Ek can be quite 
different from the SM value even under the /i — 67 constraint because of large enhancements 
of the — mixing matrix elements for the scalar operators in Eq. This correction 
of Ek leads to a change of the allowed region of the parameter 03 and eventually affects the 
possible region of AuibJ Attibj^ and A^p^B J/ipKs). 



The constraint from B(/i — 67) is somewhat model dependent. The above results depend on our choice 
of the structure of Mjv and/or T4ins- If we change these assumptions and suppress B(^ ^ e 7), the SUSY 
contribution to Am^^ can be more significant. For example, if we take the small mixing MSW solution 
for the solar neutrino anomaly, a 50% enhancement of Atob, is possible p4|. 
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In the U(2) model, the SUSY contribution to Ek can be large in a similar manner. In 
addition, there are 0(1) corrections to Mi2(-Bd) and Mi2(-Bs) so that the ^uibJ ^rriB^ and 
A^^{B — s> J/ipKs) can be different from the SM values with the same CKM matrix. 

Deviations of Ek, ^f^^B^i ^f^Bsi ^^id (pM from the SM values are plotted as functions 
of the gluino mass for tan/? = 20 in FIG. |^-^. In these figures, we fix the parameters in 
the CKM matrix as |Kife/V^b| = 0.09 and 03 = 65° and do not impose the experimental 
constraints from ek, Ams^, Am^^, and Aq^{B J/ipKs). The above features can be 
seen quantitatively in these figures. We see that ek can be different from the SM prediction 
by a factor ~ 2.5 in the SU(5) SUSY GUT with right-handed neutrinos, and the deviation is 
even larger in the U(2) model. In the mSUGRA, the deviation is smaller than 10%. As for 
Am^^ and Am^^, 0(1) deviations are possible in the U(2) model, while the deviations are 
small for the other cases. The value of (pM can also be different from the SM value (= 20i) 
significantly in the case of the U(2) model. For the mSUGRA and the SU(5) SUSY GUT 
with right-handed neutrinos, (j)M = 20i is a good approximation. 

Next, let us vary |\4fe/Vcf)| and 03 and impose all the experimental constraints explained 
above. In FIG. ^ we show possible values of A^p^B J/ipKs), AtjibJ Attib^, and 03 for 
the parameter sets satisfying the constraints. The corresponding allowed region for the SM 
is also given in each plot. Solid curves show the correlations among three quantities in the 
SM for |\4f,/V^f,| = 0.08, 0.09, and 0.10. We see that the allowed region in the SM is mainly 
determined by and Ek- 

In the mSUGRA, the deviation from the SM is not so significant since the SUSY contri- 
butions to all the Mi2's are small. 

In the SU(5) SUSY GUT with right-handed neutrinos, we see that all the allowed points 
lie between the lines corresponding to the SM values with |\4fe/Vcfe| = 0.08 and 0.10. This 
pattern arises because only the K^—K^ mixing receives SUSY corrections of 0(1) corrections 
from the SUSY loops, whereas the SUSY contributions to Mi2(-Bd) and Mi2(-Bs) are small. 
As a result, the allowed region of 03 can be shifted, and AttibJ AniB^ and A^p^B 
J/ipKs) can be different from the SM region. In other words, observables from B physics 
namely |Vu6/V^fe|, AuibJ AniB^, A^p{B J/ipKs), and 03 consistently determine a set of 
parameters in the CKM matrix in the SM analysis, though the experimental value Ek may 
not be consistent with the SM value calculated by the CKM parameters from B physics. 
The upper limit of AtubJ Aitlb^ in the plot is determined by the lower bound of Am^^ 
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FIG. 2: Ratio of ek to the SM value as a function of the gluino mass for a fixed set of the 
parameters in the CKM matrix (see text). 



21 



1.10 



1.08 - 



pa 



1.06 



pq 1.04 
1.02 



1.00 



1.08 




pq 

a" 
< 



1.06 - 



pq 1.04 



1.02 - 



1.00 



2.5 



pq 



2.0 



1.5 



pq 

S 1.0 



0.5 
0.0 



mSUGRA 
tan p = 20 



SU(5) ® 
tan p = 20 



U(2) FS 
tan P = 20 




mSUGRA 
tan P = 20 



SU(5)ev,^ 
tan p = 20 




U(2) FS 
tan p = 20 



1.10 



1.08 



1.06 



I 



3 



1.04 td 



1.02 



1.00 



1.08 



1.06 



1.04 X 



1.02' 



1.00 



2.5 




500 1000 1500 2000 

m(g) [GeV] 



500 1000 1500 2000 2500 
m(g) [GeV] 



FIG. 3: Deviations of AniB^ and Attlb^ from the SM values as functions of the gluino mass with 
the same parameter set as Fig. 0. 
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FIG. 4: Deviation of 0a/ from the SM value (= 2(j)i) as a function of the gluino mass for the U(2) 
model with the same parameter set as Fig. |2[ 



given in Eq. (|36D . 

In the U(2) model, we see that the allowed points exist outside of the region between 
iKtfc/Kbl = 0.09±0.01 lines. SUSY corrections to Mi2(-B(j) and Mi2(-Bs) in this model, unlike 
those in the SUSY GUT, are considerably large and not proportional to the corresponding 
combinations of CKM elements in the SM. Since all of Ek, AmB^/Amsj and Aqp{B — >■ 
J/ip Ks) can be corrected, there might be a mismatch in the determination of the unitarity 
triangle by the SM analysis with these quantities and 03 . In particular, we may be able to 
extract new physics contributions from observables in B physics. 

Finally, we discuss future prospects of new physics search in B decays. We expect that 
Aqp{B J/ipKs) and AtubJ Aitlb^ will be precisely measured in a few years at the B 
factories and Tevatron experiments. If we assume the SM, CKM parameters, especially 03 
can be determined from these observables with small hadronic uncertainties. By comparing 
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AJ5jf ( B ^. J/\|/ Kg ) ^3 [degree] 

FIG. 5: Scatter plots in the planes {Aqp{B J/^jjKs), AmB,,/AmBj) and (03, AmB^/ ArriBj^) 
for three SUSY models. Solid curves show the SM values with fixed |V„6/V"c6| = 0.08, 0.09 and 
0.10. 
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FIG. 6: Typical example regions. (a)Amij7AmB^ = 35 x (1 ± 0.05) and ^ J/ipKs) = 

0.75 ± 0.02. {h)AmBjAmB^ = 55 x (1 ± 0.05) and A'^^iB J/ipKs) = 0.60 ± 0.02. 
{c)AmBjAmB^ = 55 x (1 ± 0.05) and A'^^iB J/ipKs) = 0.75 ± 0.02. 

this 03 value with that derived from CP asymmetries in various B decays, we can carry out 
a consistency check of the SM and examine existence of SUSY effects. As an illustration, we 
pick up the calculated data points which satisfy the following values of A^^{B J/ipKs) 
and AniBj^niB/. 

(a) AniBjAmB^ = 35 x (1 ± 0.05), ^ J/^Ks) = 0.75 ± 0.02, 

(b) AniBjAmB, = 55 x (1 ± 0.05), ^ J/i^Ks) = 0.60 ± 0.02, 

(c) AmBjAniB^ = 55 x (1 ± 0.05), A^'^{B J/ipKs) = 0.75 ± 0.02. 

(a) corresponds to the case in which Ek, ^cpi^ ^ J/4' ^s), and AtubJ Aitlb^ are consistent 
with the SM. (b) and (c) are cases in which there are some inconsistencies among the three 
observables within the SM. The three regions are shown in FIG. ^. 
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FIG. 7: Possible region of (ps as a function of the gluino mass. 

We present the possible region of 03 in each case for the three models in FIG. |^. For 
the case (a), 03 is 60°-65° if we assume the SM. The possible value of 03 is the same as 
the SM in the mSUGRA and the SU(5) SUSY GUT with right-handed neutrinos. The 
parameter region with nig ^ 1.5 TeV in the SU(5) SUSY GUT is excluded due to the 
fi —>■ ej constraint. On the other hand, in the U(2) model, 03 can be different from the SM 
value by ~ 30° for rrig ^ 1 TeV region. 

For the cases (b) and (c), the mSUGRA (as well as the SM) is excluded because of the 
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mismatch among ek-, A^^{B J/ip Ks), and AtubJ ArriBj- In the other two models, the 
experimental value of sk can be reproduced with SUSY contributions. In the SU(5) SUSY 
GUT with right-handed neutrinos, 03 is the same as that derived from Am^^/Am^^and 
Aqp[B — s> J/ipKs) by the SM analysis. In the U(2) model, 03 can be different from the 
value of the SM. 

We have also studied the case of tan /3 = 5 and drawn the figures corresponding to FIG. |^ 
and 0. We have found that the allowed regions in these figures are almost same as those for 
tan /5 = 20. 

V. CONCLUSIONS 

In order to distinguish SUSY models by measurements at B factories, we have studied 
SUSY contributions to the - K^, Bd- Bd, and 5, - Bs mixings in three SUSY models, 
namely the mSUGRA, the SU(5) SUSY GUT with right-handed neutrinos, and the U(2) 
model. 

First, we have considered the observables Am^^, Am^^, Aqp{B J/ipKs), and Ek- 
In the mSUGRA, the deviations from the SM values are at most 10 percents for these 
observables. In the SU(5) SUSY GUT with right-handed neutrinos, the SUSY contributions 
to Ek can be large whereas those to M^i^Bd) and M^iBs) are less than 10%. In the U(2) 
model, the deviations from the SM values for all the above observables can be very large. 
In the latter two models, we may be able to see SUSY effects from the consistency check of 
the unitarity triangle among ek, Ams^/Ams^, and A^p{B J/ipKs). 

Second, we have considered cases in which the two observables /^itlbJ AmB^ and 
A™^{B — * J/ipKs) are precisely determined at the B factories and Tevatron experiments. 
We have studied how we can distinguish these different models if we determine 03 in addi- 
tion to the above two observables. We can carry out the consistency check of the unitarity 
triangle among the observables in B physics, namely AmBj^rriB^, and Aqp{B —>■ J/tp Ks), 
and 03. For the U(2) model, in particular, a large deviation from the SM value is possible. 
It is therefore very important to determine 03 precisely in theoretically clean ways from the 
decay modes, such as B ^ nn, pn, D^*'^ K^*\ D^*'^7r, D* p. These measurements are possible 
in future e~^e~ super B factories and hadron machines such as LHC-B and B-TeV. 

In this paper, we have mainly considered the consistency test of the unitarity triangle 
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through Bd decays, but there are other possibihties of finding SUSY effects in B physics. 
A new phase in the Bg — Bg mixing amphtude may affect CP asymmetries in Bg decays 
such as the Bg J/'4"P decay. These asymmetries can be measured in B experiments at 
hadron machines. For the U(2) model, these CP asymmetries could be different from the 
SM prediction Another possibility is to measure branching ratios and CP asymmetries 
in rare decays such ash sl^l~ and h — » svV. 

In conclusion we have shown that SUSY models with different flavor structures can be 
distinguished in B physics. As we have illustrated with three speciflc models, the patterns of 
the deviations from the SM in the B physics depend on the SUSY breaking mechanism and 
interactions at a high energy scale. Present and future experiments in B physics at e^e~ 
B factories and hadron machines are therefore very important to explore flavor structure of 
SUSY breakings. 
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APPENDIX A: THE CKM MATRIX AND THE UNITARITY TRIANGLE 



In this paper, we use the "standard" parameterization p2| for the CKM matrix with three 
mixing angles 9i2, 023, 9i3 and a complex phase 613. 



CKM 



v 



C12C13 



S12C13 



-S12C23 - Ci2S23Sl3e*^" C12C23 - Si2S23Sl3e*'^i=' S23C13 



S12S23 - Ci2C23Si3e 



-C12S23 - Sl2C23Sl3e 



i<5i3 



C23C13 



(Al) 



where Cij = cos 6ij and Sij = sin 6ij. The angles 0i, 02 and 03 in the unitarity triangle are 
defined as 



1 = arg 

2 = arg 

3 = arg 



In the convention ([Al|) , these angles are written in a good approximation as 



(A2a) 
(A2b) 
(A2c) 



-argMi^f(S,), 



(A3) 
(A4) 
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